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The Mars Orbiter Laser Altimeter (MOLA) detected clouds associated with 
topographic features during the polar night on Mars. While uplift generated from 
flow over mountains initiates clouds on both Earth and Mars, we suggest that the 
Martian clouds differ greatly from terrestrial mountain wave clouds. Terrestrial 
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wave clouds are generally compact features with sharp edges due to the relatively 
small particles in them. However, we find that the large mass of condensible 
carbon dioxide on Mars leads to clouds with snow tails that may extend many 
kilometers down wind from the mountain and even reach the surface. Both the 
observations and the simulations suggest substantial carbon dioxide snow 
precipitation in association with the underlying topography. This precipitation 
deposits CO2, dust and water ice to the polar caps, and may lead to propagating 
geologic features in the Martian polar regions. 

The Mars Orbiter Laser Altimeter (MOLA) used a 1.064 pm laser to measure 
distances from the spacecraft to the surface but it frequently detected reflections from 
clouds rather than the surface. In the North pole region up to 70-80% of all laser echoes 
on a given orbit have a cloud signature 1 . In some cases these echoes have distinct 
sloping edges (see Fig. 1). About half of the clouds are correlated with some sort of 
surface topography 1,2 . It has been suggested that, analogous with terrestrial mountain 
clouds, these Martian clouds contain particles too small to precipitate out and represent 
the sloping sides of propagating buoyancy or gravity waves 1,2 . Using a time dependent 
microphysical model that incorporates CO 2 microphysics 3,4 , we simulate the formation 
of CO 2 mountain wave clouds in a Martian polar atmosphere and compare the model 
results with the observations of MOLA. We find the Martian clouds contain large 
precipitating particles. 

Mountain Waves 


Orographic or mountain waves result from the flow of stratified (stable) air over surface 
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topography 5 . The upward motion of the air over the topography results in adiabatic 
cooling and thus forces buoyancy oscillations. Clouds may form within cooler regions 
of upward displacement. The sloping echo structures seen by MOLA could represent 
the tops of CO 2 clouds forming in propagating gravity waves with coherent wave fronts 
extending up from the surface (see Fig. 2). In this scenario perturbations in temperature 
induced by the vertically propagating wave lead to the formation of clouds at all 
altitudes wherever supersaturation is achieved 2 . However, at lower altitudes the 
amplitude of the propagating wave is relatively small and the resulting supersaturation 
is not sufficient to nucleate C0 2 cloud particles 6 . An alternate explanation for the 
sloping echo structures is shown in Fig. 2. The structures are still the product of 
mountain wave cloud formation. However, rather than clouds forming from the surface 
to higher altitudes along coherent wave fronts, the formation of the cloud occurs only at 
high altitude where the maximum cooling occurs 7 ' 8 . Because of the high 
supersaturations present during nucleation, condensation occurs rapidly 4,6 . Due to the 
large amount of C0 2 condensing, very large CO 2 ice particles form 4,6 . As these ice 
particles fall they are blown in the wind forming long "snow" tails. Because the air is 
nearly saturated at all altitudes and the particles are so large 4 ' 6 , even in warming 
descending air, these particles don't evaporate entirely and can fall to the ground. 

Wave Cloud Simulation 

We simulated orographic waves using a microphysical cloud model that had 
been used previously to predict cloud microphysical characterisitcs for terrestrial wave 
clouds 3 ' 4 ’ 6,7 and an analytic orographic wave model 8 ' 9 . The initial temperature profile 
was determined from the vapor pressure curve for CO 2 such that for all altitudes below 
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30 km the C0 2 humidity relative to C0 2 ice was 98%, constant with radio occultation 
measurements 10 and GCM simulations 11 . In this simulation small (r 0 ~ 1 pm) ice 
covered dust grains are assumed to be the nucleation sites for the C0 2 cloud 
particles 4,6,12 . 

The surface topography is modulated with time giving the appearance of a 
parcel of air moving at a fixed speed over the surface. The size and frequency change 
of the topography along with the surface wind speed will determine the properties of the 
wave induced 5,8 . In general, for stable stratified conditions, wide ridges (~ 10 km) with 
wind speeds of less than 5 m sec' 1 provide the most favorable conditions for the 
formation of vertically propagating waves. Under these conditions the greatest amount 
of cooling occurs at levels just below the wave breaking altitude 9 and a large number of 
new C0 2 cloud particles form. 

For cloud particles formed in the simulation, the back-scattering amplitude, 
described by the scattering phase function, is solved for using Mie theory. From this 
back-scatter amplitude the back-scatter cross section, ( 3 , of the simulated clouds can be 
calculated. The back-scatter cross section is defined here as the product of the average 
particle cross section with the back-scattering amplitude. Typical back-scatter cross 
sections are 100 - 300 pm 2 km" 1 . To compare these back-scatter cross sections with 
MOLA observations, the observed average back-scatter cross section was calculated 
using the single scatter approximation 13 . The calculated back-scatter cross sections for 
a MOLA pass over the northern polar region are shown in Fig. lb. In any single cloud 
event, the particle back-scatter cross sections are fairly constant, with an average cross 
section of 200 pm 2 km' 1 corresponding to a particle radius of approximately 1000 pm 
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(assuming 35% of the atmospheric mass at 20 km is condensed to form a cloud with an 
optical depth over a 400 meters path length equal to one at a wavelength of 1 .064 pm). 
However, from one pass to another, the cross section may vary by as much as a factor of 
two. 

Comparison with MOLA observations 

The simulation presented here uses the measured topography from the MOLA 
pass number 226. This topography is indicated in Figure 3 by the gray shaded regions. 
For a 5 m sec' 1 surface wind flowing from the left to the right a cloud forms over the 
first ridge at an altitude of about 1 8 km. Before reaching the ground the falling CO 2 ice 
crystals blow approximately 20 km down wind. The back-scatter cross section 
calculated from the simulated cloud are plotted (contours) together in Fig 3 along with 
the MOLA measured echo returns (filled dots). There is good agreement with both the 
altitude and the horizontal extent of the clouds, as well as the location of formation 
relative to the surface topography. The color shading of the MOLA echo data indicates 
the magnitude of the back-scatter cross sections calculated from the MOLA 
observations for the same portion of pass 226. Here too there is good agreement 
between the magnitude of the back-scatter cross section, as calculated from the cloud 
simulation, and the back-scatter cross sections derived from MOLA observations. 

In many of the simulations, including the one just presented, the cloud particles 
fall to the ground before evaporating. The calculated integrated mass flux of CO 2 ice 
particles to the surface for a simulation of MOLA pass number 266 is shown in Fig. 3 
(green dashed line belonging to the right vertical axis). The average size of the ice 
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particles that are being deposited to the surface is about 200 pm. The ice coated dust 
particles that these CO 2 ice crystals formed on are typically about 5 pm in diameter 
giving an integrated CO 2 dust to ice mass ratio of < 0.01 %. This is about 5 times lower 
than estimates using the Thermal Emission Spectrometer observations to model the 
radiative behavior of the polar caps 14 . The discrepancy occurs because of further 
modification of the ice crystals once on the surface 15 and by a seasonal enhancement of 
dust concentrations as CO 2 ice evaporates in the summer months. 

Conclusions 

The formation of carbon dioxide clouds requires a supersaturation of about 
35% 6 . This requirement largely determines the characteristics of CO 2 cloud particles 
and where they will form 3 ' 4 . The adiabatic cooling associated with orographic waves is 
sufficient to achieve the supersaturations required to nucleate new cloud particles. The 
wave-breaking altitude primarily determines the altitude of the clouds since that is 
where the largest amplitude waves and greatest temperature perturbations occur. 

Significant amounts of CO 2 ice can be deposited to the surface during one of 
these snowstorms. In the simulations done here over an one-hour period as much as 
0.75 grams of CO 2 snow (Fig 3) may be deposited per cm 2 . The deposition of snow will 
be associated with the location of the wave clouds and hence the underlying 
topography. Snow accumulation at the tops and lee sides of ridges may explain some of 
the differences in surface brightness measured by TES and seen by the Mars Orbital 
Camera (MOC). In particular, bright crater rims, plateaus, and linear streaks running 
the length of ridges 14, 16 may be the result of enhanced accumulation of CO 2 and water 
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ices associated with the wave activity produced by these topographic features. The 
continued deposition of dust and retention of water ice associated with the nuclei of the 
CO 2 cloud particles will further modify the polar terrain by burying seasonally 
condensed C0 2 (down wind of the wave). Over geologic time, this process may result 
in a propagation of terrain, analogous to dune fields and may assist in the modulation of 
terraces and layering in specific regions. Being optically thick, these carbon dioxide 
clouds most certainly effect the radiation field during the polar night and lead to low 
brightness temperatures as suggested by Forget and Pollack and others ’ ’ ’ . 

While mountain wave clouds are a common occurrence on Earth, the high 
condensate mass on Mars make these clouds unlike their terrestrial counterparts. 
Terrestrial wave clouds are composed of relatively small particles (around 5 pm in 
size), which do not fall far and evaporate quickly 7 . The snow generated in Martian 
wave clouds is able to survive for a long time as it is blown downwind. The long "tail" 
of snow originating at altitudes in excess of 5 km can extend for 10's of kilometers 
down wind and deposit substantial amounts of ice to the surface. These clouds 
represent a subset of clouds within the polar hood. Diffuse, optically thin water clouds 
composed of fine particles (< 2 pm) are ubiquitous at all Martian latitudes 3, 22, 23 and 
additional CO 2 clouds can form if ice crystals are lifted by synoptic scale systems into 
slightly supersaturated regions. Regardless, because of their frequency, total mass and 
spatial extent, these carbon dioxide wave clouds play a critical role in both the CO 2 
mass and radiation budgets of the poles. 
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Methods 

Microphysical Cloud Model 

The microphysical model includes the processes of nucleation, 
condensation/evaporation, and sedimentation. It is based on terrestrial models that have 
been used to study cirrus wave clouds, polar stratospheric clouds, and stratus clouds 1 . A 
contact parameter of 0.95 was used for all simulations. Transport of gases, particles and 
potential temperature from one altitude bin to another is allowed via sedimentation, 
diffusion and vertical advection. Sedimentation velocities are calculated for two 
regimes. A Stokes-Cunningham expression is employed for Reynolds numbers < 10' 2 
and empirical interpolations are used for Reynolds numbers greater than 10‘ 2 . Diffusion 
is handled by an eddy diffusion scheme in which the diffusion coefficients are 
calculated by a boundary layer turbulent mixing model. Advection is calculated using a 
linear wave model. 

The model has 50 altitude bins each 1 .0 km in thickness. The size distributions 
of both ice-coated dust and cloud particles are divided across 50 bins that range in size 
from .01 pm to 5000 pm. The ice coated dust distribution is assumed to initially have a 
log normal size distribution with a modal radius of r 0 = 0.1 pm (note on top of pg. 4 you 
say ro is 1 pm)and an integrated number concentration of 50 cm' 3 mixed uniformly with 
height. 


Wave Model 

For a horizontal wind speed, constant in the vertical and in time, the equation for 
the vertical velocity perturbation, w', associated with internal gravity waves is 
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In this equation x and z are the horizontal and vertical distances respectively and N is 
the Brunt-Vaisala frequency given by 

f(r,-r) V 

where F/ and Fare the dry and atmospheric lapse rates respectively. Typical values for 
N in the Martian atmosphere are between 5 x 10' 3 and 1 0‘ 2 s' 1 . A solution to this 
equation for statically stable stationary waves forming over a ridge is 

x 

w(x, z ) = e !/2H Re(w T exp[/(^x + m s zj]) 


yielding the dispersion relation m 2 = N 1 ju 1 - k 2 . Here k and m are the horizontal and 
vertical wave numbers respectively and subscript s indicates the s‘ h Fourier component. 
The atmospheric scale height is H = RT/Mg . For mechanical topographic forcing the 
gravity wave vertical velocity amplitude w s can be related to the height of the 
topography (at the zonal wave number k s ) hr by: 

w, = k s uh T 

A vertical propagating gravity wave will grow in amplitude as p Aa (where p is pressure) 
until the waves become unstable and "break" 7 . The breaking altitude of these waves 
occurs at the level where the local lapse rate, the sum of that due to both the waves and 
the mean state, first exceeds the dry adiabatic lapse rate 8 . Above the breaking altitude 
the wave dissipates through turbulence. 


Back-scatter Cross Sections 

The back-scatter cross section is a measure of the ability for a particle to scatter light 
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backward and depends on the size, or cross section, of the scattering particle and its 
shape. For a single scatter approximation the back-scatter cross section, /?, may be 
written 


/3 =4n 


P, A, T 


where P, and P r are the laser power transmitted and received respectively. A, is the area 
of transmitted laser footprint at the cloud tops, A r is the collecting area of the receiver, / 
is the range distance to the cloud tops, and T is the attenuation of the laser by the 
atmosphere. The atmospheric attenuation is defined as T = ae~ 2r where a is the surface 
albedo and T is the atmospheric total optical depth. 
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Figure 1 Echo returns from MOLA pass number 232 (a) and the average back-scatter 
cross section calculated from laser energy return at 1 .067 pm for the same MOLA pass 
(b). The long sloping features that extend from the surface to nearly 6 km in altitude 
(part a) are the subject of the discussion here and are believed to be carbon dioxide 
"snow tails". 

Figure 2. A cartoon showing the stream lines that would be produced from a stationary, 
altitude independent flow over a single ridge. Also shown in this figure are the regions 
of possible cloud formation (indicated by the shaded regions). As the waves propagate 
upward their amplitude will increases leading to greater cooling and thus higher 
supersaturations (5) at higher altitudes. It is in these regions that CO 2 clouds will form. 
Precipitation from these clouds will fall downward, blowing down wind (indicated by 
the speckled region). (This fig looks odd, it seems to have a big black block for the 
snow. Maybe it is just not read right by my computer.) 

Figure 3. The simulated average back-scatter cross section at 1.067 pm for a 5 m sec' 1 
surface flow over the topography as measured by MOLA during pass number 226. The 
shaded region represents the measured topography. MOLA echoes for the same pass 
are shown as solid colored dots. The color of each dot represents the back-scatter cross 
section (color scale at the top of the figure) derived from MOLA laser energy data as 
described in the text. Contours show the back-scatter cross sections calculated for the 
simulation. The contour magnitude is designated by the same color as the MOLA 
echoes (dots). Also shown (green dashed line with right axis) is the calculated mass 
flux of CO 2 ice particles to the surface for a segment of MOLA pass 226 between 0 and 
100 km. For the simulation presented here in an one hour snow storm approximately 
0.75 g cm' 2 of CO 2 snow may be deposited to the surface. 
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